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Abstract 

A  thermochemical  kinetics  analysis  has  been  carried  out  for  growth  on  the 
(100)-(2x  1):H  diamond  surface  using  a  pair  of  previously  proposed  mechanisms  which 
operate  sequentially.  Half  of  the  growth  on  such  a  surface  is  accounted  for  by  inser¬ 
tion  into  dimer  bonds,  while  the  other  half  is  accounted  for  by  addition  across  troughs 
between  dimer  bonds.  The  latter  mechanism  is  slower  and  therefore  controls  the  over¬ 
all  growth  rate  on  this  surface.  This  result  can  explain  the  success  that  the  latter 
mechanism  has  had  in  predicting  growth  rates  in  a  variety  of  systems.  We  suggest 
that  growth  at  step  sites  is  favored  on  steric  and  thermochemical  grounds  and  can 
account  for  atomically  smooth  surfaces  on  diamond. 
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Introduction 


In  the  past  few  years  several  detailed  growth  mechanisms  have  been  proposed  to 
explain  the  chemical  vapor  deposition  (CVD)  of  diamond  films  [1-10].  In  general, 
these  models  have  used  the  formal  resemblance  between  the  bonding  and  structure 
in  diamond  and  in  alkanes  to  postulate  [6]  that  chemistry  on  the  diamond  surface 
could  be  understood  in  terms  of  the  very  well  known  chemistry  of  alkanes.  In  effect, 
the  assumption  is  made  that  the  chemistry  of  a  diamond  surface  is  controlled  by  the 
local  electronic  environment,  as  is  true  for  alkanes,  and  that  bulk  band  structure  and 
surface  states  play  no  role.  If  this  postulate  is  valid,  then  considering  all  that  is  known 
about  alkane  chemistry,  we  may  be  able  to  understand  CVD  diamond  growth  at  a 
level  of  detail  beyond  that  for  any  other  CVD  process. 

Considerable  research  in  the  field  has  been  aimed  at  determining  the  “growth” 
species,  i.e.,  the  gas  phase  species  directly  responsible  for  diamond  formation  [11, 
12],  and  recent  experiments  have  demonstrated  that  the  CHz  radical  is  the  primary 
growth  species  in  the  CVD  systems  that  have  been  examined  [13-19].  One  proposed 
chemical  kinetics  mechanism  [2],  which  takes  C//3  as  the  growth  species,  successfully 
predicts  [20-22],  experimental  growth  rates  for  both  RF  [21]  and  DC  [23]  plasma 
torches,  for  flames  at  low  and  atmospheric  pressure  [15, 24,  25],  and  for  filament  sys¬ 
tems  as  a  function  of  pressure  and  composition  [16, 17],  without  the  use  of  adjustable 
parameters.  Its  predictions  are  compared  to  experiment  in  Figure  1.  The  ability 
of  the  model  to  predict  the  relative  growth  rates  for  these  very  different  systems  is 
striking.  However,  we  have  pointed  out  [2]  that  the  near-perfect  absolute  agreement 
is  fortuitous — our  estimated  uncertainty  is  a  couple  of  orders  of  magnitude  in  either 
direction  [2].  This  is  because  there  is  uncertainty  not  only  in  gas  phase  rate  and 
thermodynamic  parameters  themselves,  but  also  in  how  these  parameters  should  be 
converted  for  use  with  surface  kinetics.  As  discussed  previously  [2],  the  mechanism 
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models  the  diamond  surface  with  the  9-carbon  molecule  bicyclononane  (BCN),  which 
in  effect  takes  the  diamond  surface  to  be  an  unreconstructed  (100)  surface.  However, 
it  is  unlikely  that  such  a  surface  would  be  stable  because  of  very  large  H  ~  H  steric 
repulsions  [26].  Recent  STM  and  AFM  experiments  have  shown  [27-29]  that  diamond 
surfaces  can  be  rather  rough  on  an  atomic  scale,  but  at  least  some  of  the  surface  re¬ 
constructs  to  the  (100)-(2x  1):H  form  [28,30,31].  On  this  surface  each  carbon  atom 
is  bonded  twice  to  carbons  in  the  bulk,  once  to  a  surface  hydrogen,  and  once  to  a 
neighboring  surface  carbon  making  a  “dimer”  and  forming  a  5-membered  ring,  as 
seen  in  Figure  2.  A  quantitative  analysis  of  the  steric  forces  on  various  (100)  surfaces 
shows  that  this  dimer  reconstruction  substantially  reduces  surface  stress  compared 
to  the  unreconstructed  surface  [26,32].  In  this  paper  we  propose  that  the  growth 
kinetics  on  the  (100)-(2x  1):H  surface  can  be  understood  by  combining  two  differ¬ 
ent  mechanisms — a  “trough”  mechanism  (identical  to  the  BCN  mechanism)  and  a 
“dimer”  mechanism — which  operate  sequentially  on  the  two  types  of  sites  present  on 
this  surface.  This  picture  of  growth  can  explain  the  success  that  the  BCN  mechanism 
hais  enjoyed. 


Analysis 


Any  acceptable  diamond  growth  mechanism  must  be  both  fast  enough  to  explain 
observed  growth  rates  and  thermodynamically  favorable  enough  so  that  the  overall 
reaction  in  the  reverse  direction  (etching)  is  negligible  [2,33,34].  The  latter  require¬ 
ment  comes  from  the  experimental  observation  that  etching  by  atomic  hydrogen  is 
extremely  slow.  In  this  section  we  describe  our  analysis  for  the  thermochemistry  and 
kinetics  of  proposed  growth  mechanisms  which  satisfy  these  conditions. 


3 


Structure  of  the  Model  Compound 

The  first  step  is  to  choose  a  model  structure  to  represent  the  diamond  crystal  and 
surface.  Use  of  very  large  model  structures  is  a  requirement  if  steric  energies — which 
can  be  quite  large — are  to  be  calculated  accurately.  This  is  because  the  calculation 
must  take  into  account  the  fact  that  lattice  atoms  close  to  a  surface  reaction  site  may 
move  somewhat  to  relieve  steric  repulsions,  but  it  must  also  take  into  account  that 
the  crystal  as  a  whole  is  extremely  rigid.  For  estimating  steric  energies,  we  used  the 
crystal  slab  shown  in  Figure  2,  which  is  1.6  by  1.6  by  0.6  nm  (8  layers)  thick  and 
contains  330  carbon  atoms  and  20  surface  dimer  bonds.  Reaction  takes  place  near 
the  center  of  the  slab.  Although  this  crystal  is  relatively  large,  when  finding  the  steric 
energy  it  was  necessary  to  hold  twenty  atoms  at  the  base  and  edges  of  the  crystal 
fixed  in  order  to  prevent  the  crystal  from  flexing  due  to  the  large  surface  tensile  stress 
induced  by  the  dimers.  For  estimating  heat  capacities  and  entropies  we  used  a  smaller 
(160  carbon  atoms)  l.O  by  0.6  by  0.6  nm  thick  diamond  slab  because  of  limitations 
of  the  computer  program. 

Gas  Phase  Species  Concentrations 

For  the  models  discussed  here  we  assume  that  growth  occurs  from  reaction  of 
the  CH3  radical  at  the  diamond  surface.  We  take  gas  phase  concentrations  from 
measurements  and  modeling  [12,  35-37].  For  our  typical  filament-assisted  growth 
conditions  with  the  substrate  temperature  1200  K  and  the  pressure  20  torr  we  have 
Xh  =  2  X  10~^  and  Xchs  =  2  x  lO'**,  where  Xi  is  the  mole  fraction  of  species  i. 

Rate  Constants 

Rate  constants  for  the  gas-surface  reactions  that  appear  in  these  mechanisms  have 
not  been  measured.  In  order  to  estimate  their  values  we  have  made  the  assumption 
[4]  that  the  reaction  cross  section  per  surface  site  is  equal  to  the  reaction  cross  section 
per  equivalent  site  in  a  prototype  gcis  phase  reaction.  According  to  our  analysis  for 
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an  abstraction  by  a  gas  phase  species  of  mass  nig  this  assumption  leads  to  [4] 

ks  =  kgU^lfi/nigY^^, 

where  fc,  and  kg  are  the  surface  and  gas  phase  rate  constants,  respectively,  is  the 
reduced  mass  of  the  reactants  in  the  prototype  gets  phase  reaction,  and  the  symmetry 
number  n,  is  the  number  of  identical  sites  on  the  surface  structure  being  examined. 
In  other  words,  the  surface  rate  constant  equals  the  gas  phase  rate  constant  corrected 
for  the  effect  of  collision  frequency  and  for  symmetry.  We  use  the  same  approach 
for  radical  recombination  reactions,  with  one  difference.  We  consider  the  radical  site 
on  the  gas  phase  molecule  to  be  equivalent  to  two  surface  radical  sites  because  the 
nearly-planar  gas  phase  radical  has  two  equally  probable  sides  for  reaction,  whereas 
the  surface  radical  has  only  1  side  from  which  reaction  can  occur.  Thus,  for  radical 
recombination  we  have 

ks  =  O.bnskglfi/nigY^^. 

Assuming  a  constant  reaction  cross  section  neglects  some  factors  which  may  affect 
the  surface  rate  constant.  Steric  hindrances  and  reaction  barriers  might  differ  on  the 
surface  from  those  of  the  prototype  gas  phase  reactions.  Also,  transition  state  theory 
shows  that  the  pre-exponential  factor  is  determined  by  the  entropy  change  A5^  in 
forming  the  activated  complex;  even  if  the  reaction  potential  surface  is  the  same  for 
the  prototype  and  surface  reactions,  differences  in  the  translational  and  rotational 
contributions  to  may  alter  the  pre-exponential  factor  somewhat.  These  effects 
are  ignored  here,  since  they  are  difficult  to  estimate  accurately  and  are  not  expected 
to  qualitatively  affect  our  conclusions. 

The  rate  constants  are  shown  in  Tables  1  and  2.  (Definitions  of  the  species  are 
given  below.)  Two  surface  isomerization  reactions  in  Table  1,  G  and  //,  are  treated 
differently.  We  do  not  expect  a  significant  barrier  for  Reactions  G,  so  we  have  cissigned 
it  a  large  rate  constant  which  simply  insures  that  this  reaction  is  in  partial  equilibrium. 
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The  rate  constant  for  Reaction  H  has  been  estimated  by  Musgrave,  Goddard,  and 
Harris  [38]  using  ah  initio  quantum  techniques. 

Thermochemistry 

We  estimate  values  for  the  enthalpy  change  A//  and  entropy  change  A5  in  order 
to  determine  AG  =  A//  —  TAS'  for  each  reaction.  From  AG  we  can  calculate  rate 
constants  for  reverse  reactions,  which  are  related  to  the  forward  rate  constants  by 

/ie,  =  kj/kr  = 

where  A'e,  is  the  equilibrium  constant.  (Thermochemical  quantities  are  tabulated 
for  standard  states  of  1  atmosphere,  whereas  we  are  using  concentration  units.  For 
reactions  in  which  the  number  of  moles  changes,  kr  must  be  multiplied  by  a  factor 
which  converts  between  atmospheres  and  moles/cm^.) 

To  estimate  thermodynamic  parameters  we  have  suggested  the  use  of  molecular 
mechanics  [34].  With  this  technique,  a  force  field  is  established  that  describes  the 
interaction  of  every  atom  in  the  system  with  every  other  atom.  In  general,  both 
the  functional  forms  for  the  interactions  and  the  associated  coefficients  are  chosen 
empirically.  A  number  of  force  fields  have  been  proposed  for  carbon  atoms  [39-41],  but 
the  most  widely  used  and  extensively  tested  force  fields  for  obtaining  thermodynamic 
quantities  are  MM2  [42]  and  MM3  [43],  which  we  use  in  this  paper.  These  force 
fields,  which  were  designed  specifically  to  estimate  steric  repulsions  in  highly  strained 
alkanes,  give  each  bond  a  set  of  force  constants  (e.g.,  stretching,  bending,  torsion) 
for  displacement  from  “natural”  lengths  and  angles,  and  give  each  atom  van  der 
Waals  attractions  and  repulsions  to  the  other  atoms  in  the  system.  Our  codes  use 
the  MM2  or  MM3  force  fields  to  adjust  the  position  of  every  atom  in  the  system  in 
order  to  minimize  the  total  strain  energy  which  is  the  result  of  bonds  bending, 

stretching,  and  twisting  in  response  to  van  der  Wcials  attractions  and  repulsions.  The 
heat  of  formation  Hj  for  a  molecule  is  obtained  by  combining  this  calculated  Entrain 
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with  bond  enthalpies  Hbond  calculated  from  an  associated  group  additivity  scheme. 
MM2  and  MM3  predict  H j  of  even  very  crowded  and  highly  strained  molecules  to 
within  typically  0.5  kcal/mole  (0.02  eV/molecule),  which  is  not  much  greater  than 
the  experimental  error  [42,43].  On  the  other  hand,  since  MM3  is  a  purely  empirical 
force  field  and  has  not  been  calibrated  against  a  diamond  surface,  we  do  not  expect 
the  predicted  energetics  to  be  so  accurate  for  the  present  work.  For  example,  we 
estimate  that  fixing  atoms  at  the  base  and  edges  of  the  model  compound  (see  above) 
introduces  an  uncertainty  of  around  1  or  2  kcal/mole  in  relative  heats  of  formation. 
(Only  relative — not  absolute — heats  of  formation  appear  in  the  analysis.)  The  MM3 
force  field  can  also  calculate  vibrational  frequencies.  For  stable  structures  we  use  the 
enthalpy  calculations  from  MM3.  Since  the  radical  force  field  in  that  code  has  not  yet 
been  finalized  [44],  we  use  the  MM2  force  field  to  estimate  the  difference  in  the  heats 
of  formation  of  stable  species  and  radical  species  with  a  hydrogen  removed.  Heats 
of  formation  calculated  by  MM2  and  MM3  differ  significantly  only  where  there  are 
non-bonded  H-H  distances  less  than  0.2  nm  [43],  which  occur  for  the  species  C^M 
and  CqHM.  For  these  cases,  an  error  of  1  to  2  kcal/mole  may  be  introduced  by  our 
use  of  MM2.  The  entropy  is  calculated  primarily  by  using  the  MM3  code.  However, 
we  have  also  included  a  symmetry  term  — /Jinn,  to  take  into  account  the  presence 
of  identical  atoms  for  some  of  the  surface  structures,  and  we  treat  -CHj  and  -CH^ 
groups  as  free  rotors.  Treating  these  as  free  rotors  is  a  good  approximation,  since  the 
barriers  to  rotation  are  found  to  be  typically  less  than  1  kcal/mole.  Thermodynamic 
quantities  calculated  at  1200  K  for  the  two  mechanisms  considered,  relative  to  the  C5 
starting  structure  (see  below),  are  shown  in  Table  3  and  4. 

Reaction  Mechanisms 

Carbon  atoms  can  add  onto  the  (100)-(2xl):H  surface  at  two  types  of  sites,  in¬ 
serting  into  dimer  bonds  or  adding  across  troughs  between  dimer  bonds.  To  complete 
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a  new  monolayer,  half  of  the  carbons  must  add  to  the  dimer  sites  and  half  to  the 
trough  sites.  The  added  carbons  may  then  dimerize  with  their  neighbors,  forming  a 
new  (2x1)  surface  on  a  new  layer,  with  dimer  rows  orthogonal  to  the  original  rows. 
Recently,  Garrison  and  co-workers  [45]  have  proposed  a  mechanism  for  insertion  into 
dimer  bonds  which  involves  exclusively  reactions  that  have  low  activation  barriers  and 
well-known  analogs  in  hydrocarbon  chemistry.  (Another  proposed  mechanism  for  in¬ 
sertion  into  dimer  bonds  [46]  includes  a  step  with  a  3-center  transition  state.  Since  the 
barrier  at  that  transition  state  is  calculated  [46]  to  be  very  high,  it  is  not  considered 
further.)  For  addition  across  the  troughs  we  use  the  bicyclononane  (BCN)  mecha¬ 
nism,  which  Weis  originally  proposed  [2]  to  explain  growth  on  the  unreconstructed 
(100)  surface.  Although  the  atoms  making  up  a  trough  are  isomorphic  to  those  on 
an  unreconstructed  (100)  surface,  on  the  (100)-(2xl):H  surface  the  strain  energy  as¬ 
sociated  with  this  structure  is  substantially  reduced.  Thus,  the  BCN  mechanism  can 
be  applied  more  realistically  to  troughs  on  the  (100)-(2x  1);H  surface. 

Dimer  mechanism 

The  steps  in  the  dimer  mechanism  of  Garrison  et  al.  [45]  are  reviewed  in  Table  1. 
This  mechanism  begins  with  the  dimer  structure  on  the  (100)-(2x  1):H  surface,  which 
we  denote  Cs  (Figure  3a)  because  of  the  presence  of  5-membered  rings.  In  the  initial 
step  (Reaction  A)  a  surface  hydrogen  is  abstracted  by  a  gas  phase  H  atom  forming 
the  tertiary  surface  radical  Cl  and  the  //|“*  molecule.  Either  (Reaction  B)  or 
C (Reaction  C)  can  recombine  at  this  radical  site  to  give  C5  or  C5M ,  respectively. 
A  hydrogen  atom  from  the  methyl  group  of  C5M  can  be  abstracted  to  give  C5M' 
(Reaction  D),  which  in  turn  may  react  with  (Reactions  E  or  F).  In  the  next 
step  (Reaction  G)  C5M*  isomerizes  to  CJ  (Figure  3b),  a  species  which  contains  both 
a  radical  site  and  a  double  bond.  Reaction  G  is  called  a  0  scission  reaction  [47] 
because  the  C-C  bond  two  carbons  away  (or  “/?”)  from  the  radical  site  is  broken 
during  this  isomerization  reaction.  (Two  other  C-C  bonds  can  also  break  via  a  0 
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scission  reaction;  products  of  those  reactions  are  not  considered  here.)  The  radical 
carbon  in  Q  can  react  with  either  end  of  the  carbon-carbon  double  bond.  Addition  to 
the  end  closest  to  the  bulk  diamond  reverses  Reaction  6",  while  addition  to  the  other 
end  forms  the  species  (Reaction  H,  Figure  3c).  Finally,  Cg  and  Ce  interconvert 
by  Reactions  /  and  J.  We  note  that  a  total  of  two  hydrogen  abstraction  reactions 
are  necessary  to  form  Ce  from  Ce- 

Trough  Mechanism 

The  steps  involved  in  adding  carbon  across  the  trough  on  the  (100)-(2x  1):H 
surface  are  identical  to  those  proposed  previously  in  the  trough  mechanism  [2,  3] 
(except  that  formation  of  long-chain  hydrocarbons,  Reaction  j,  is  not  considered 
here),  and  they  are  shown  in  Table  2.  Briefly,  a  methyl  radical  adds  at  a  radical  site 
formed  from  abstraction  of  one  or  both  of, the  trough  hydrogens  (see  Figure  3e).  This 
step  is  followed  either  by  abstraction  of  the  other  trough  hydrogen  and  one  of  the 
methyl  hydrogens  (Reactions  d  or  f  and  k  after  Reaction  c)  or  by  abstraction  of  a 
methyl  hydrogen  (Reaction  /  after  Reaction  s)  to  give  the  CeB  structure,  which  has 
a  carbon  atom  that  bridges  the  trough.  For  reasons  discussed  below,  we  allow  the 
trough  mechanism  to  occur  only  adjacent  to  Ce-type  sites  just  formed  by  the  dimer 
mechanism,  as  shown  in  Figure  3e.  Thus,  the  starting  structure  is  called  CeBHs 
which  refers  to  the  two  H  atoms  pointing  toward  each  other  across  the  trough  and 
next  to  a  newly-formed  Ce  species.  The  final  bridged  structure  is  called  CeB.  We 
note  that  a  total  of  three  hydrogen  abstraction  reactions  are  necessary. 

Results 


Dimer  mechanism 

The  time  dependence  for  the  dimer  mechanism  is  calculated  by  integrating  the 
rate  equations  for  Reactions  A  through  J.  Initially,  the  fraction  of  the  dimer  sites 
which  are  Cs  is  1.0,  and  the  fraction  of  all  other  species  in  the  mechanism  is  0.  The 
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calculation  shows  the  surface  rapidly  reaches  steady  state,  within  less  than  100  ms. 
However,  there  are  no  irreversible  steps  in  this  mechanism,  since  those  steps  which 
have  a  large  enthalpy  reduction  also  have  a  substantial  reduction  in  entropy.  The 
result  is  that  the  reaction  sequence  is  reversible,  in  conflict  with  experiment. 

This  situation  can  be  easily  remedied  by  adding  Reactions  /\ ,  L,  and  M  to  the 
mechanism,  which  creates  a  dimer  on  a  new  layer,  which  we  call  Layer  n  +  1.  The 
resulting  structure,  C&HH ,  is  shown  in  Figure  3e.  C^HH  is  generated  when  a 
species  is  formed  directly  in  front  of  or  behind  (eclipsing)  a  C^T*  species,  which  has 
a  radical  site  at  the  top  carbon  (Figure  3d), 

Ce  +  CeT*  +  =  CeHH  +  ■  (M) 

Because  AG  0  for  Reaction  M,  it  is  effectively  irreversible  (assuming  that 
we  have  not  omitted  some  other  etching  reactions  of  C&H H  which  could  convert  it 
back  to  Cs),  and  C^H H  ultimately  covers  half  a  monolayer.  The  rate  of  formation 
of  that  half-monolayer  would  correspond  to  an  effective  linear  growth  rate  of  around 
7  microns/hour.  However,  since  after  approximately  100  ms  growth  is  effectively 
stopped  until  the  trough  sites  can  be  filled  in,  the  actual  growth  rate  is  limited  by 
whichever  process  is  slower. 

Trough  Mechanism 

Although  the  dimer  mechanism  including  Reactions  K  through  M  is  relatively 
fast  and  irreversible,  it  accounts  for  growth  only  on  dimer  sites;  the  trough  mech¬ 
anism  adds  carbon  atoms  at  the  sites  which  separate  the  dimers.  We  take  as  the 
starting  point  for  the  trough  mechanism  the  structure  C^HH  produced  by  the  dimer 
mechanism,  that  is,  one  dimer  already  formed  on  Layer  n  -f  1.  Because  of  the  large 
negative  values  of  AG  for  Reactions  k  and  /,  formation  of  C^B  is  irreversible  via  the 
trough  mechanism.  Furthermore,  CsB  cannot  react  back  to  C^HH  via  the  reverse  of 
the  dimer  mechanism  since  formation  of  a  dimer-type  bond  across  the  trough  is  not 
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energetically  feasible.  C&B  formation  is  followed  by  a  rapid  reactions  analogous  to  Re¬ 
actions  K  through  M  of  the  dimer  mechanism  in  which  a  pair  of  eclipsed  CeB  species 
reacts  to  form  a  dimer  bond  on  Layer  7i  +  1.  Thus,  the  combined  dimer/trough  mech¬ 
anism  starts  with  two  dimer  bonds  and  creates  a  new  surface  containing  two  dimer 
bonds,  one  each  from  the  dimer  and  trough  parts  of  the  mechanism. 

There  are  a  couple  of  approximations  implicit  in  this  picture  that  should  be 
pointed  out.  First,  in  order  to  keep  the  reaction  mechanism  manageable  we  have 
made  an  independent  sub-unit  approximation,  assuming  that  the  kinetic  and  ther¬ 
modynamic  parameters  for  each  dimer-trough  pair  are  independent  of  the  status  at 
neighboring  dimer-trough  pal  j.  Second,  we  note  that  when  all  of  the  dimer-trough 
pairs  have  reacted  the  new  surface  is  again  a  (100)-(2x  1);H  surface,  but  it  has  not 
been  returned  to  a  structure  identical  to  the  starting  structure  because  the  dimer 
bonds  formed  on  Layer  n  -f  1  are  at  right  angles  to  those  on  Layer  n.  In  fact,  because 
of  the  abcabc . . .  structure  of  the  diamond  lattice.  Layer  n  is  not  reproduced  until 
Layer  n  -f  4.  However,  the  structures  of  each  of  these  layers  are  congruent — they 
are  merely  shifted  or  rotated.  Thus,  we  assume  that  the  growth  rate  of  each  layer 
between  n  and  n  +  4  takes  place  at  the  same  rate. 

The  most  important  difference  between  the  two  mechanisms  from  the  point  of 
view  of  kinetics  is  that  the  trough  mechanism  is  considerably  slower,  due  partly  to 
the  additional  H  atom  abstraction  that  must  occur  with  the  trough  mechanism.  The 
steady  state  lineai  growth  rate  of  the  combined  mechanism  is  calculated  to  be  0.5 
/im/hr,  limited  by  the  trough  portion  of  the  mechanism.  This  rate  is  controlled  by 
two  factors.  The  first  factor  is  the  rate  of  addition  of  CHj  to  the  surface,  which  is 
determined  by  the  fraction  of  the  surface  sites  which  are  radicals  and  by  the  rate 
constants  for  Reactions  c  and  s.  The  second  factor  is  the  fraction  of  the  CH3  radicals 
which  desorb  before  they  can  be  incorporated  into  the  crystal.  For  example,  the 
reverse  of  Reaction  c,  in  which  the  C H3  group  from  C%H M  desorbs,  is  several  times 
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faster  than  Reactions  d  or  /,  in  whicii  the  surface  growth  processes  continue.  In 
contrast,  the  reverse  of  Reaction  s,  in  which  the  C H3  from  Ce  M  desorbs,  is  an 
order  of  magnitude  slower  than  the  irreversible  formation  of  C^B  in  Reaction  1.  This 
difference  is  a  consequence  of  the  considerably  greater  steric  repulsion  in  CeH  M 
compared  to  C^M ,  which  leads  to  spontaneous  C desorption  lifetimes  at  1200  K 
of  100  //s  and  11  ms  for  the  two  species,  respectively.  As  a  result,  most  of  the  CqB 
which  is  formed  passes  through  Reaction  s  rather  than  Reaction  c,  even  though  Cg* 
is  a  relatively  high  energy  species.  The  hydrogen  atom  concentration  also  plays  a 
role,  since  irreversible  incorporation  of  carbon  into  the  lattice  requires  abstractions 
by  gas  phase  H  atoms  which  compete  with  methyl  desorption.  We  have  published 
previously  our  predictions  for  the  effect  of  H  atom  concentration  on  growth  rates  [3, 
4],  but  we  are  not  aware  of  any  experimental  data  which  test  them. 

Discussion 


Since  a  (100)-(2x  1):H  surface  contains  equal  numbers  of  dimer  and  trough  sites, 
the  dimer  and  trough  mechanisms  ultimately  contribute  equally  to  growth.  However, 
steric  repulsion  controls  how  and  when  these  mechanisms  play  their  roles.  The  inter¬ 
actions  between  these  mechanisms  can  be  discussed  with  reference  to  the  drawing  in 
Figure  4,  which  labels  dimer  sites  as  a,  c,  e,  and  and  trough  sites  ^ls  6,  d,  /,  and  h. 

Growth  of  Layer  n  -|-  1  can  commence  only  with  the  dimer  mechanism  and  not 
the  trough  mechanism.  This  is  because  spanning  a  trough  giving  a  hypothetical  C5B 
structure  would  create  highly  strained  bonds  and  would  be  unstable  because  the  C-C 
bonds  would  have  to  be  too  long.  Thus,  the  first  reactions  are  reversible  formation 
of  Ce  species  on  the  dimer  sites  a,  c,  e,  or  g  followed  by  the  irreversible  formation  of 
CqHH.  CqH H  sites  will  not  in  general  exist  simultaneously  at  a  and  c  or  at  e  and 
g  because  that  would  lead  to  large  H-H  steric  repulsions  across  the  trough.  This 
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would  shift  the  reactions 


CsM-  ^  c:  ^  c; 

strongly  to  the  left  and  slow  the  dimer  mechanism.  Therefore,  once  a  new  dimer 
bond  CeH H  forms  at  c  and  for  example,  reaction  at  a  and  e  to  make  another 
CeH H  structure  there  is  unlikely.  However,  the  presence  of  a  C&H H  species  at  c  and 
g  shortens  the  distance  across  the  trough,  which  allows  troughs  ft,  d,  /,  and  h  to  be 
spanned.  Thus,  C^B  structures  are  formed  with  the  trough  mechanism,  and  these 
react  with  each  other  to  form  additional  dimer  bonds  on  Layer  +  1.  At  this  point, 
with  the  trough  hydrogens  eliminated,  the  dimer  mechanism  rapidly  adds  carbon  at 
a  and  e.  Thus,  growth  on  the  (100)-(2x  1):H  surface  occurs  in  a  coordinated  fashion, 
first  up  and  down  rows  of  (eclipsed)  dimers,  then  up  and  down  adjacent  rows  of 
(eclipsed)  troughs,  then  up  and  down  the  remaining  rows  of  (eclipsed)  dimers. 

According  to  our  analysis,  the  steady  state  growth  rate  at  dimer  sites  is  con¬ 
siderably  faster  than  at  trough  sites.  Since  a  complete  new  layer  cannot  form  until 
addition  has  taken  place  at  each  trough  site,  the  trough  mechanism  is  rate  limiting. 
This  result  may  provide  an  explanation  for  the  observation  that  experimental  growth 
rates  can  be  calculated  with  the  trough  mechanism,  as  seen  in  Figure  1. 

Atomically  Smooth  Surfaces  on  Diamond 


In  this  work  we  have  represented  the  surface  by  an  ensemble  of  independent  sub¬ 
units  consisting  of  a  pair  of  dimer  bonds  next  to  a  pair  of  troughs,  such  as  sites  a,  ft, 
e,  and  /  (see  Figure  3a).  To  some  extent  this  representation  is  adequate.  For  exam¬ 
ple,  the  difference  in  steric  energy  between  Cb  and  C^M  changes  by  only  around  1 
kcal/mole,  depending  on  whether  or  not  a  C^HH  structure  is  adjacent  to  the  Cs.  Fur¬ 
thermore,  the  independent  sub-unit  representation  does  allow  us  to  predict  a  certain 
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amount  of  coordination  between  subunits.  For  example,  because  of  the  alternating 
nature  of  the  dimer/trough  mechanism,  we  found  that  neighboring  dimer  rows  cannot 
grow  until  the  trough  between  them  has  been  spanned.  However  this  approach  does 
not  account  for  the  relatively  smooth  surfaces  that  have  sometimes  been  observed  over 
long  ranges  on  the  (100)-(2x  1):H  surface  [28,30,31].  On  metals  and  some  other  types 
of  crystals,  smooth  surfaces  are  the  result  of  loosely  bound  adatoms  diffusing  to  steps 
and  kinks,  which  lowers  their  energy  (by  increasing  their  coordination  numbers)  and 
leads  preferentially  to  growth  at  those  sites.  This  process  produces  fast  lateral  growth 
and  smooth  surfaces.  In  contrast,  it  is  unlikely  that  chemically  bound  hydrocarbon 
adsorbu,tes  would  be  mobile  on  a  hydrogenated  diamond  surface  [48].  Furthermore, 
the  coordination  number  of  hydrocarbon  adsorbates  is  unaffected  by  whether  they  are 
bound  at  step  or  terrace  sites.  Thus,  diffusion  of  adsorbates  to  steps  and  kinks  would 
not  provide  a  mechanism  for  explaining  atomically  smooth  growth  on  diamond. 

We  suggest  instead  that  steric  and  thermochemical  factors  can  favor  growth  at 
step  sites  over  growth  on  smooth  surfaces.  For  example,  we  have  already  shown 
that  growth  on  atomically  smooth  (111)  surfaces  is  highly  unlikely  with  a  straight¬ 
forward  abstraction-addition  mechanism  [34].  As  we  have  pointed  out  above,  smooth 
unreconstructed  (100)  surfaces  may  well  not  exist,  but  we  found  previously  [3]  that 
methyl  radical  addition  at  step  sites  with  (100)  character  can  readily  occur  because 
these  sites  are  in  general  less  sterically  crowded.  Similarly,  addition  to  (110)-  or 
dimer/trough-type  sites  may  be  sterically  easier  at  steps  than  on  flat  surfaces.  If 
growth  intermediates  are  more  stable  at  step  and  kink  sites,  growth  would  occur 
preferentially  there  for  the  sai  .3  reason  that  growth  with  the  trough  mechanism  oc¬ 
curs  preferentially  through  C^M  rather  than  through  CsHM.  The  result  would  then 
be  fast  lateral  growth  leading  to  smooth  surfaces. 
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Table  1 


DIMER  MECHANISM 


REACTION 

^Cj200 

it* 

A. 

Cs  + 

Cs*  +  i/|" 

-14.2 

2.52  X  10^4^-7^00/ RT 

B. 

Cs*  + 

Cs 

-59.5 

1.0  X  10^3 

C. 

Cs*  +  Ci/f  * 

CsM 

-30.2 

5.0  X  10^2 

D. 

CsM  + 

CsM- +  //f" 

-16.9 

2.81  X 

E. 

CsAT  + 

CsM 

-57.6 

1.0  X  10^3 

F. 

CsAT  + 

Cs*  +  C//|“ 

-27.4 

3.0  X  10’3 

G. 

CsM- 

Q 

+4.3 

1.0  X  10*3 

H. 

Q 

Cs* 

-21.2 

6.9  X  ioi2e-8ioo//?T 

I. 

Cs*  + 

Cs 

-48.5 

1.0  X  10*3 

J. 

Cs  + 

Cs*  +  7/f " 

-25.2 

1.26  X  10*^e-73Oo/«T 

A'." 

Cs  + 

^  CeT  +  " 

-22.9 

9.0  X 

L.‘= 

CeT  + 

^  Cs 

-50.8 

1.0  X  10*3 

M.^ 

CsT*  +  Cs  + 

^  CsHH  +  Hf* 

-61.1 

1.8  X  lO^T^e- 3000/^^7' 

a.  units  are  kcal/mole 

b.  units  are  cm^,  moles,  seconds 

c.  Analogous  reactions  can  be  written  for  Cs- 
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Table  2 


TROUGH  MECHANISM 


REACTION 

AC»2“ 

k 

a.C&HH 

Ceir  +  ^1“ 

-16.1 

1.3  X  lOi^e-^aoo/flr 

b.  Ceir  + 

CeHH 

-57.6 

1.0  X  10'3 

c.  Ce^  +  C^|“" 

CeHM 

-20.6 

5.0  X  10^2 

d.CeHM  + 

C^M  +  //|“* 

-30.3 

1.3  X  ioJ4e-73oo//?T 

e.  C^M  + 

CeHM 

-64.0 

1.0  X  10*3 

/.  CeHM  + 

Ce/ZAT  +  //f* 

-19.1 

2.8  X  io7T2e-77°o/«^ 

5.  Ceff  AT  + 

CeHM 

-54.6 

1.0  X  10*3 

/i.  Ce^Af* 

Ce’M 

-11.2 

1.0  X  10* 

i.  CsHM"  + 

CeH'  +  C^|“ 

-34.0 

3.0  X  10*3 

ib.  Ce/^AT  + 

Cefi  +  ^1“*  ' 

-67.3 

1.3  X  I0*4e-7300//ir 

/.  Ce'M  + 

— 

CeB  +  //f  * 

-56.1 

2.8  X 

m. 

Ce*^  +  //f* 

-22.3 

1.3  X  i0*4e-^3oo/fi7’ 

n.  Q/f  + 

CeHH 

-51.4 

1.0  X  10*3 

0.(76//*  + 

;:e 

Ce"  + 

-18.7 

1.3  X  l0*'‘e-^3oo/Hr 

p.  Ce**  +  //®“" 

CeH- 

-55.0 

1.0  X  10*3 

9.  Ce'H  + 

Ce"  + 

-12.5 

4.5  X  io6r2e-50oo/«r 

r.  Ce**  +  //3“* 

Ce^H 

-61.2 

1.0  X  10*3 

5.  Ce**  +  C^f* 

Ce-M 

-32.2 

5.0  X  10*2 

t.  Ce//* 

Ce^H 

-6.2 

1.0  X  10® 
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Table  3 


THERMOCHEMISTRY  FOR  DIMER  MECHANISM 


SPECIES 

kcal/mole 

5“ 

cal/mole-Kelvin 

kcal/mole 

fls 

Cs 

0.0 

0.0 

0.0 

2 

Cs* 

+44.0 

-0.1 

+44.1 

1 

CsM 

9.1 

+20.5 

-15.4 

3 

CsAT 

+51.3 

+20.4 

+26.8 

2 

Q 

+56.2 

+20.9 

+31.1 

1 

Ce* 

+27.9 

+15.0 

+9.9 

1 

Ce 

-5.8 

+  14.5 

-23.2 

2 

CeT 

+32.8 

+  17.2 

+  12.2 

1 

56.6 

34.3 

+  15.4 

6.4 

41.1 

-42.9 

46.1 

62.8 

-29.3 

a.  Values  for  H  and  S  and  G  are  for  1200  K.  Values  for  all  surface  species  referenced 
to  C$. 
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Table  4 


THERMOCHEMISTRY  FOR  TROUGH  MECHANISM 


SPECIES 

//“ 

kcal/mole 

5“ 

cal /mole- Kelvin 

kcal/mole 

ra* 

CeHH 

+  11.4 

+21.0 

-13.8 

1 

CeH‘ 

+51.7 

+  19.5 

+28.4 

1 

Ce‘H 

+46.1 

+  19.9 

+22.2 

1 

cr 

+90.0 

+  18.4 

+68.0 

1 

CeHM 

+26.9 

+40.3 

-21.5 

3 

CeHAT 

+65.8 

+40.1 

+  17.7 

1 

+53.3 

+38.9 

+6.5 

3 

CeB 

+40.9 

+26.9 

+8.7 

1 

a.  Values  for  H  and  S  and  G  are  for  1200  K.  Values  for  all  surface  species  referenced 
to  C5. 
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Captions 


1.  Comparison  of  measured  growth  rates  with  growth  rates  predicted  from  nu¬ 
merical  simulations  [20-22]  using  the  trough  mechanism  [2].  The  environments 
simulated  were  an  oxygen-acetylene  flame  [24]  (solid  squares),  a  220  torr  DC 
arcjet  [23]  (open  circle),  an  atmospheric-pressure  RF  torch  [21]  (open  squares), 
a  40  Torr  acetylene-oxygen  flat  flame  [25]  (solid  triangle),  and  a  hot-filament 
system  [22]  (solid  circle).  The  calculations  use  the  parameters  given  in  the  pa¬ 
per  proposing  the  BCN  mechanism  [2]  with  the  exception  of  the  rate  constant 
for  H  atom  recombination  at  radical  sites.  This  value  was  adjusted  from  10*'*  to 
5  X  10^^  cm^/mole-s  because  the  former  value  implies  a  recombination  rate  equal 
to  twice  the  collision  rate. 

2.  The  diamond  slab  used  as  the  model  compound  for  calculating  steric  repulsion 
energies. 

3.  Surface  structures  considered  in  this  work.  The  letters  >4,  B,  C,  and  D  represent 
adducts  on  the  crystal  slab. 

(a)  Cs  {A  =  H);  Q  {A  =');  C^M  {A  =  C//3);  C,M‘  {A  = 

(b)  Q. 

(c)  Ce*  [B  =•);  Ce  {B  =  H). 

(d)  Ce  eclipsing  a  CeT*. 

(e)  C^HH  {C  =  H,D  =  H);  {C  =  H,D  =*); 

CeHM  {C  =  H,  D  =  CHs);  QM  (C  =*,  D  =  CH^)- 
CeHM*  (C  =  i/,  D  =  CH;)-,  Cg**  (C  =*,  D  =*); 

Ce//*//  (C  =*,  D  =  H). 

(f)  CeB. 
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4.  Structure  of  a  portion  of  the  (100)-(2x  1):H  surface.  Positions  labeled  a,  c,  e, 
and  g  are  sites  for  addition  with  the  dimer  mechanism  positions  labeled  6,  d,  /, 
and  h  are  sites  for  addition  with  the  trough  mechanism. 
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